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Claim of Priority 

This application claims priority from U.S. Provisional Patent 
Application Serial Number 60/465,975 bearing Attorney Docket Number 
12003007 and filed on April 28, 2003. 

Field of the Invention 

This invention relates to using external polymeric analytical techniques to 
predict in-vivo polymeric performance. 

Background of the Invention 

Polymers have transformed material science in the past seventy years. 
Medical science has benefited from the rheology and morphology of polymers 
for the construction of prosthetic medical devices and implants. 

Any medical device for mammalian usage requires extensive screening of 
performance to assure safety and efficacy of use. Stages of clinical trials 
proceeds from animal patients to human patients. The opportunity cost of such 
clinical trials is significant in terms of risk, expense, and time. 

Summary of the Invention 

The art needs a means to limit the effort of clinical trials for polymer 
performance in medical devices. 

The present invention solves that problem by providing a method of 
using external polymeric analytical techniques to predict in-vivo polymeric 
performance. 
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"External" means beyond the mammalian body. "In-vivo" means within 
the mammalian body. 

"Polymeric" means of or relating to "polymers", which can be naturally- 
occurring polymers, synthetically-produced polymers, combinations of 
5 naturally-occurring polymers and synthetically-produced polymers by physical 
or chemical interaction, or combinations thereof 

The ability to predict performance is extremely important in the 
development of in-vivo biomedical devices, such as implants for skeletal joints. 
Often the devices contain a polymeric component that provides damping as well 
10 as structural integrity. The devices are exposed to complex deformations (both 
shear and compression-extension) and complex loading schedules (cyclic loads 
of varying amplitude/magnitude). 

Aside from bio-compatibility considerations, the lifetime of the 
polymeric component is typically determined by two characteristics of the 
1 5 polymeric component. 

The first is accumulated, irreversible strains (micro-failures). 
Accumulated micro-failure is usually characterized by fatigue testing (lifetime 
vs. load or deformation amplitude). 

The second "failure" mode involves the viscoelastic deformational 
20 response of the polymeric component to the loading schedule. 

The present invention focuses on this second "failure" mode. 
Specifically, the present invention uses time-temperature superposition to define 
the stress- and time-dependent properties of a polymeric component, extending 
to 20 years of in-vivo usage, from predictions based on external laboratory 
25 measurements requiring only a few weeks. Boltzmann's superposition is used 
to define the deformation response to a cyclic loading schedule. Lifetime is 
assessed by the ability of the device to maintain a required geometry during the 
imposed loading schedule. 
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One aspect of the present invention is a method for predicting in- vivo 
polymer performance in a mammalian body via external experimentation, 
comprising the steps of: 

(a) determining stress levels and the duration times of stress application, 
5 denominated a daily stress cycle that the polymeric component will be subjected 

to during mammalian body activities for the lifetime of the device; 

(b) determining compressive creep and recoverable creep compliance at 
selected stresses, as determined in step (a), and temperatures to obtain master 
curves of compressive creep compliance that analytically define compressive 

1 0 creep strain as functions of time and stress; 

(c) determining single creep/creep recovery strain curves for each 
mammalian body activity of step (a) calculated from the compressive creep 
compliance master curves of step (b), wherein a creep phase is interrelated to a 
recovery phase, to obtain a compressive creep model for a single stress cycle; 

1 5 (d) determining a total residual strain after repeated stress cycles from 

the compressive creep model of step (c), to obtain a compressive creep model 
for the polymer subjected to repeated daily mammalian activities, denominated 
compressive stress cycles; 

(e) computing maximum compressive strain curves from a compressive 
20 strain during an activity in relation to number of days of activity, wherein the 

number of days can range from a few days to several years; and 

(f) using the maximum compressive strain curves of step (e) to predict 
in- vivo performance of the polymer by calculating polymer deflection as a 
function of original geometry of the device to be made of the polymer. 

25 A feature of the present invention is that the polymer need not be placed 

in- vivo without some estimation of its performance in- vivo. 

An advantage of the present invention is that the external prediction can 
reduce patient risk, duration of clinical trial product screening activities, and 
attendant expenses associated with both. 
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The present invention can achieve the following goals: 

1 . Demonstrated Utility of Viscoelastic Property Characterization 
in Performance Modeling of Polymeric Biomaterials. 

2. Time-Temperature Superposition Can Often Be Employed to 
Provide Long-Time Properties (Temperature-Acceleration). 

3. Boltzmann's Superposition Provides Basis for Defining 
Stress/Strain Response for Various Loading Histories/Requirements. 

4. Performance Modeling Can Provide Basis for Material Selection. 
Other features and advantages will be revealed when amplifying on the 

embodiments of the invention in comparison with the following drawings. 

Brief Description of the Drawings 

Fig. 1 shows a typical Fatigue S-N Failure Curve. 

Fig. 2 shows a depiction of a sample geometry subjected to load. 

Fig. 3 shows a graph comparing stress and strain over time, showing 
creep and creep recovery. 

Fig. 4 displays the creep compliance for various compressive stresses 
(loads) for the polymeric component. 

Fig 5 displays creep compliance measurements that were obtained at a 
constant stress of 145 pounds per square inch (psi) at test temperatures of 37, 
50, 70, 85 and 100°C. 

Fig. 6 shows a resulting master curve for the 37°C creep compliance. 

Fig. 7 depicts the strain response for a single stress cycle. 

Fig. 8 is a graph of a compressive creep model for a single stress cycle 
of 158 psi. 

Fig. 9 shows the accumulation of residual strain as a function of number 
of days or number of stress cycles. 

Fig. 10 shows the actual external strain response for each activity for the 
first two days of the device's predicted in-vivo lifetime. 
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Fig. 1 1 shows the predicted accumulated "residual" strain for each 
activity, calculated during a "resting" period (end of day), as a function of 
number of days into the predicted in-vivo lifetime of the disc device. 

Fig. 12 shows the time-dependent strains of the method of the present 
5 invention. 



Embodiments of the Invention 

10 In- Vivo Polymer Predictive Performance 

Polymers have found considerable popularity in biomedical applications 
because of their unique properties. Polymers possess both "solid-like" and 
"liquid-like" behavior. They can be used as "carriers" that change from a liquid 
to solid, depending on their environment. They can provide structural 

15 properties (rigidity) as well as energy dissipation (mechanical damping). They 
can provide friction reduction at moving surfaces. 

Some of the biomedical applications that involve polymers include drug- 
delivery systems, dental adhesives, tissue adhesives, support fabric, prosthetic 
devices, skeletal implants, etc. 

20 Each of these non-limiting in-vivo examples requires that the polymer 

provides the designed performance properties while being relatively inert to its 
surrounding, i.e. bio-compatible. 

The present invention uses viscoelastic property characterization to 
define the capability of a polymer to provide the "designed performance 

25 properties". (Bio-compatibility is not addressed.) 

For example, one item benefiting from the present invention is a spinal 
implant device that replaces a damaged disc. The disc device contains a 
polymer that provides both structural integrity and mechanical damping. The 
in-vivo device is exposed to complex deformations involving both shear and 

30 compression-extension in response to a complex loading schedule, i.e. cyclic 



5 



File No. 1200307R 

METHOD OF PREDICTION OF IN- VIVO POLYMER PERFORMANCE BY EXTERNAL 
EXPERIMENTATION 

stresses that vary in magnitude and duration. The major "designed performance 
property" is that the device has to provide structural integrity in terms of extent 
of deformations for various loads for a lifetime of 30 years. 

Lifetime predictions for "structural integrity" of a polymeric component 
5 are typically assessed by two techniques: cyclic fatigue and long-time creep. 
Fatigue studies emphasize failure by accumulated irreversible strains (micro- 
failures) or accumulated damage. Accelerated testing is accomplished by 
testing at cyclic stress levels that are higher than those involved in the 
application. Table 1 below summarizes this point. 

10 



Table 1 


Fatigue 


Lomg-Time Creep 


Cycles to Failure vs. Load 


Viscoelastic Property 
Characterization 


Accumulated Damage 


Model Deformation vs. Time and 
Load History 


Accelerated by Increasing Load 


Accelerated by Increasing 
Temperature 


Viscoelasticity 


Irreversible Deformations 


Survival Without Excessive Deformation 



An S-N curve is obtained that describes the lifetime as a function of 
cyclic stress level. Fig. 1 shows a typical Fatigue S-N Failure Curve. If the 
polymer component "survives" the anticipated number of cycles that the device 
1 5 incurs during a 30-year lifetime at a cyclic stress level that exceeds the 

anticipated in- vivo stress, then the device is "acceptable". Although the S-N 
curve is generally interpreted in terms of accumulated damage, the change in 
cyclic deformation is also influenced by the viscoelasticity of the polymer 
component. 
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Returning to Table 1, long-time creep studies emphasize the viscoelastic 
properties of the polymer. Properties are measured to provide the capability to 
calculate/model the deformation of the polymeric component as functions of 
time and loading history for the projected lifetime of the device. Performance 
5 "acceptability" is directly related to deformation of the polymeric component 
throughout the proposed lifetime. Typically, creep and creep recovery 
experiments are performed at various constant stress levels to determine 
"linearity" of response. Irreversibility and non-linearity can be defined but 
accumulated damage is not part of the model. Acceleration of testing is 

1 0 accomplished by increasing the temperature of testing. 

Fig. 2 shows the geometry of the polymeric component in the disc 
device that will serve as the "model" for the device and for the compressive 
creep/creep recovery test. The cylindrical sample is confined at its top and 
bottom surfaces by an adhesive. The "fixed" boundary of the polymer causes a 

1 5 very complex stress/strain distribution throughout the sample (as can be shown 
by Finite Element Analysis, a technique known to those skilled in the art) 
compared with the essentially bi-axial strain that evolves from the 'lubricated" 
boundary compressive creep test specimen. As shown in Fig. 2, a constant 
compressive stress (load divided by cross-sectional area) is applied to the 

20 sample at time zero and the time-dependent strain (Ah divided by initial height, 
H, or thickness of the sample) is measured. At some time t r the stress is 
removed from the sample and the time-dependent strain (creep recovery) is 
measured. Fig. 3 shows that time-dependent response graphically. 

Table 2 shows five equations concerning creep compliance and 

25 viscoelastic performance of the polymeric material. 



Table 2 

I D(t) = e(t)/a Creep Compliance 
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II 


D = D(t) only; "Linear" Viscoelastic 


III 


D = D(t, a) "Non-linear" Viscoelastic 


IV 


D creep — ^creep recovery R©V©TSibl© 


V 


Dcreep > D C reep recovery "Non-reversible" 



The compressive creep compliance, D(t), is defined as the ratio of the 
time-dependent strain divided by the constant stress. If D(t) is independent of 
stress, then the behavior of the polymer is classified as "linear" viscoelastic. If 
5 D(t) depends upon stress, then the behavior is classified as "non-linear" 
viscoelastic. The creep recovery portion of the experiment allows the 
calculation of a recoverable creep compliance, D r (t). If D r (t) is equal to D(t), 
then the deformation is considered to be "reversible". If Dr(t) is less than D(t), 
then the deformation is considered to be "irreversible". "Irreversibility" can 

1 0 originate from a micro-fracture or from a purely viscous response. 

Fig. 4 displays the creep compliance for various compressive stresses 
(loads) for the polymeric component. As noted, the compliance is dependent on 
stress and is thus considered "non-linear". As applied stress increases the 
compliance decreases. Specifically, strain increases less than proportionally 

15 with increasing stress, i.e. the polymeric component "strain-hardens", which 
indicates the presence of a finite network structure. Each test required 
approximately 1 1.5 days at a test temperature of 37°C. Even though there is 
considerable information available within this creep compliance data, the 1 1.5 
days of time is insufficient to use for predicting behavior to 30 years (10957 

20 days or 9.5x1 0 8 s). 

Time-Temperature Superposition is often used to develop a master curve 
of properties, such as creep compliance, for extrapolating to long-times. 
Increasing the testing temperature effectively accelerates the creep process so 
that it occurs in a shorter time- frame. 
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Fig. 5 displays creep compliance measurements that were obtained at a 
constant stress of 145 psi at test temperatures of 37, 50, 70, 85 and 100°C. It is 
noteworthy that the compliance curves have regions of similar curvature and 
that compliance increases with increasing testing temperature. The acceleration 
5 of the creep process with increasing temperature allows one to shift each curve 
along the log-time scale to superpose with the data for the temperature of 
interest, i.e. for a temperature of 37°C. 

The resulting master curve for the 37°C creep compliance is shown in 
Fig. 6. The use of Time-Temperature Superposition (temperature acceleration) 
10 provided creep compliance characterization at 37°C for times in excess of 10 10 s 
(300 years) from creep measurements at five temperatures involving creep times 
of only 55 hours for each test. 

At this juncture, creep compliance of the polymeric component of the 
disc device has been determined as a function of compressive stress for times up 
15 to 30 years. This permits external prediction of in- vivo performance. 

Time-dependence is provided by the 37°C compliance master curves. 
Temperature-dependence (not necessary for the present application) is provided 
by the a T function. Stress-dependence is provided by compliance master 
curves at the stresses of interest. The completion of the performance prediction 
20 model requires a knowledge of the loading schedule that is imposed on the 
device during the in- vivo lifetime. 

Table 3 shows a typical daily load/stress schedule for the polymeric 
material in use. 



Table 3 


Activity 


Characteristic 
Stress (psi) 


Daily Duration of 
Activity (hours) 


Supine 


58 


8.4 


Sit/Walk/Stand 


158 


14 
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Table 3 


Activity 


Characteristic 
Stress (psi) 


Daily Duration of 
Activity (hours) 


Jumping/Playing 


238 


1.5 


Bending/Lifting 


558 


0.1 



Although the device could be subjected to other than compressive 
stresses, the present discussion models only the deformational response to 
5 compressive stresses. As shown, the device is continuously subjected to a 
compressive stress of at least 58 psi throughout its lifetime. Each activity has 
an associated incremental increase in stress for a characteristic duration time. 

(Certain liberties are taken to assign multiple types of activities to 
characteristic stress levels. The actual stress schedule will influence the 
10 deformational response but the trends that will be displayed by the present 
model will not be significantly affected.) 

The modeling of the deformation of the device then uses Boltzmann's 
superposition as described in the equations shown in Equation VI below. 

15 VI: e(t) = aiD(t-ti) + a 2 D(t-t 2 ) + . . . + a n D(t-t n ) 

where: <j\ = compressive stress (constant) 
tj = time of stress application 
n = number of stress level changes 

20 Simply stated, the total deformational response of a polymer is equal to 

the sum of the individual strain responses for each incremental applied stress 
from the beginning of the loading/stress history and extending to the time of 
interest. Specifically, there is are deformational response contributions from the 
first stress application (first term) and each additional stress application that 

25 persist throughout the lifetime of the device. 
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Fig. 7 depicts the strain response for a single stress cycle. Notice the 
persistence of a "residual" strain after removal of the stress. Mathematically, 
the time-dependent strain response from a single stress cycle can be represented 
by a two-term equation, as shown at the bottom of this slide. For a load/stress 
schedule having 3 daily stress cycles, the model describing time-dependent 
strain for 2 days involves 1 3 terms is shown in Table 4. 

Non-linearity of the compressive creep model is introduced by selecting 
the appropriate master curve representation for compliance that is associated 
with the magnitude of the stress for the activity. The model describing time- 
dependent strain for a 30-year lifetime involves 65,746 terms, as seen in Table 
4. 



Table 4 


Deformation Model: 4 Daily Activities; 2 Days 


s(t) = cnD,(t) 


Activity 1 (continuous) 


+ a 2 D 2 (t-t 2 ) - cr 2 D 2 (t-t 3 ) 


Activity 2/Day 1 


+ a 3 D 3 (t-U) - a 3 D 3 (t-t 5 ) 


Activity 3/Day 1 


+ a 4 D 4 (t-t 6 ) - a 4 D 4 (t-t 7 ) 


Activity 4/Day 1 


+ a 2 D 2 (t-t 8 ) - cr 2 D 2 (t-t 9 ) 


Activity 2/Day 2 


+ a 3 D 3 (t-t 10 ) - a 3 D 3 (t-tn) 


Activity 3/Day 2 


+ a 4 D 4 (t-ti 2 ) - a 4 D 4 (t-ti 3 ) 


Activity 4/Day 2 






For 30 year lifetime: 
6 x 365.25 x 30 + 1 


65,746 terms! 



Fortunately, and unexpectedly, the number of terms can be reduced to 
only 7 to provide an accurate description of the time-dependent strain. Consider 
the strain response to a single stress cycle, e.g. the Sitting/Walking/Standing 
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activity as seen in Fig. 8. The strain for a daily stress cycle at day two is the 
sum the response at day two (creep recovery of day one stress cycle) and the 
response at day one (second stress cycle at day two). The strain at day three 
would be the sum of the responses at day one and at day two and at day three. 
5 And so on. 

Fig. 9 shows the accumulation of residual strain as a function of number 
of days or number of stress cycles. Thus for each activity there is a "residual" 
or accumulated strain contribution. Strain at any time within the cycle can be 
considered as the sum of the residual strains of each activity (from previous 

1 0 stress cycles) plus the strain resulting from the current stress cycle. 

Fig. 1 0 shows the actual external strain response for each activity for the 
first two days of the device's predicted in- vivo lifetime. It should be noted that 
the activities that are associated with the higher stress levels have the shorter 
duration. Although the strain during the high stress application is high, the 

1 5 recovery process is relatively rapid after removal of the stress. Activities that 
have longer times produce larger amounts of "residual" strain. 

Fig. 1 1 shows the predicted accumulated "residual" strain for each 
activity, calculated during a "resting" period (end of day), as a function of 
number of days into the predicted in- vivo lifetime of the disc device. The total 

20 strain during resting is the sum of each "residual" strain component for each 
activity. It is noteworthy that the supine activity is associated with the lowest 
stress level, but its response contributes to nearly 76% of the total deformation 
during resting. The bending/lifting activity involves the highest stress levels, 
but because of its short duration the residual strain from this activity is a 

25 negligible component of the total strain during resting. It is also noteworthy 

that if "micro-damage" occurs in the actual application, it is likely caused by the 
higher-stress-level activity, even if duration is short. This aspect is not included 
within the present model. 
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The maximum strain that is produced within the disc during any activity 
is essentially comprised of the total residual strain for each activity (the strain 
response to previous activities, i.e. stress history) and the creep deformation 
caused by the given activity. 
5 Fig. 12 shows the time-dependent strains that are associated with each of 

the activities, satisfying the present invention. At this juncture, one can assess 
and predict external the amount of deflection that will occur in- vivo, within the 
device as a function of activity-type and the number of days into the device's 
lifetime. It has been found using the method of the present invention that the 
10 disc deflection is related to strain through the original geometry of the device. 
Based upon the requirements of the device to be used in-vivo, one can assess 
via external analysis whether the polymer component provides acceptable 
predictive performance before clinical trials are initiated. 



15 Usefulness of the Invention 

The present invention provides viscoelastic property characterization for 
performance modeling of bio-medical devices that incorporate a polymeric 
component and are load-bearing during service. 

Time-Temperature Superposition can often be used to accelerate testing 
20 of pertinent properties. Boltzmann's superposition provides a mathematical 
methodology for determining the time-dependent strain that develops in 
response to an imposed stress history. The modeling of the present invention 
provides an opportunity to describe and predict behavior of the device during 
in-vivo service as well as it providing a basis for evaluating alternate 
25 "candidate" polymers for use in the construction of the device. 

The invention is not limited to the above embodiments. The claims 

follow. 
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